Variable TeV emission as a manifestation of jet formation in M87? 
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ABSTRACT 

It is proposed that the variable TeV emission observed in M87 may be pro- 
duced in a starved magnet ospheric region, above which the outflow associated 
with the VLBA jet is established. It is shown that annihilation of MeV pho- 
tons emitted by the radiative inefficient flow in the vicinity of the black hole, 
can lead to injection of seed charges on open magnetic field lines, with a den- 
sity that depends sensitively on accretion rate, n± oc m 4 . For an accretion rate 
that corresponds to the inferred jet power, and to a fit of the observed SED by 
an ADAF model, the density of injected pairs is found to be smaller than the 
Goldreich- Julian density by a factor of a few. It is also shown that inverse Comp- 
ton scattering of ambient photons by electrons (positrons) accelerating in the gap 
can lead to a large multiplicity, ~ 10 3 , while still allowing photons at energies 
of up to a few TeV to freely escape the system. The estimated gap width is not 
smaller than 0.01r s if the density of seed charges is below the Goldreich- Julian 
value. The VHE power radiated by the gap can easily account for the luminosity 
of the TeV source detected by H.E.S.S. The strong dependence of injected pair 
density on accretion rate should render the gap emission highly intermittent. We 
also discuss briefly the application of this mechanism to Sgr A*. 



1. Introduction 

Combined VLBA and TeV observations of M87 reveal a rapidly varying TeV emission 
that appears to be associated with the m.a.s VLBA jet (Acciari et al. 2009). The rapid 
flaring activity of the TeV source, with timescales t = liday day as low as 1-2 days, implies 
a source size of d ~ 4r s t day for a black hole mass M BH = 4 x 10 9 solar masses^ 



1 While early results based on M87 gas kinematics suggested a black holes mass Mbh = (3.2±0.9) x 1O 9 M0 
(Macchettto et al. 1997), recent modeling indicates M B h = (6.4±0.5) x 1O 9 M (Gebhardt & Thomas 2009). 
For the following, a fiducial value of 4 x 10 9 is adopted for M87. 
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This, and the fact that the TeV emission appears to be correlated with the VLBA jet but 
not with emission from larger scales (and in particular HST-1), motivates the consideration 
that the observed TeV photons originate from the black hole magnetosphere (Neronov & 
Aharonian 2007; Rieger & Aharonian 2008). A plausible magnetospheric process discussed 
in the literature is curvature and/or IC emission by particles, either hadrons or leptons, 
accelerating in a vacuum gap of a starved magnetosphere. An alternative explanation for 
the observed fast TeV variability is emission from small regions located at larger radii, 
r ~ 100r 9 , as, e.g., in the misaligned minijets model of Giannios et al. (2010), or perhaps 
interaction of the jet with red giant stars closer to its base (Barkov et al. 2010). 

The presence of the VLBA jet implies that a force-free (or ideal MHD) flow is established 
on scales < 100r 9 (Walker et al. 2008, Acciari et al. 2009), so that the magnetosphere is 
anticipated to be screened in the sense that the invariant E • B nearly vanishes everywhere. 
However, as in pulsar theory, there must be a plasma source that replenishes charges which 
escape the system (both, to infinity and across the horizon) along the open magnetic field 
lines in the polar region. The nature of this plasma source is poorly understood at present. 

The injection of charges into the magnetosphere may be associated with the accretion 
process. Direct feeding seems unlikely, as charged particles would have to cross magnetic 
field lines on timescale shorter than the accretion time in order to reach the polar outflow. 
Free neutrons that may be produced in a radiative inefficient accretion flow (RIAF) can 
cross field lines, however, they will decay over a distance ~ 0.03r s for a 10 9 M Q black hole, 
and even if existent at sufficient quantity, will not reach the inner regions. On the other 
hand, MeV photons that are emitted by the hot gas near the horizon can annihilate in the 
polar region to produce charged leptons. Below, it is shown that the density of the charges 
thereby injected depends sensitively on the accretion rate and the conditions in the RIAF. 
Naive estimates suggest that in case of M87 this process cannot provide complete screening 
at accretion rates that corresponds to the inferred jet power, and to a fit of the observed SED 
by an ADAF model. Those estimates are, however, highly uncertain, as explained below. 

Another plausible plasma source is cascade formation in starved magnetospheric re- 
gions. The size of the gap then depends on the conditions in the magnetosphere and the 
pair production opacity. As shown elsewhere (Levinson 2000), vacuum breakdown by back- 
reaction is unlikely, as it requires magnetic field strength in excess of a few times 10 5 G, 
higher than the equipartition value for Eddington accretion. Pair production via absorption 
of TeV photons by the ambient radiation field is more likely. However, the spectrum of the 
VHE photons observed by H.E.S.S. extends up to ~ 10 TeV (e.g., Aharonian et al. 2006), 
and the assumption that these photons originate from the magnetosphere (or even the VLBA 
jet) implies that the pair production opacity at these energies must not exceed unity. This 
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raises the question whether pair cascades in the magnetosphere can at all account for the 
multiplicity required to establish a force-free flow. 

In what follows it is shown that inverse Compton scattering (IC) of ambient photons 
by electrons (positrons) accelerating in the gap can lead to a large multiplicity, ~ 10 3 , while 
still allowing photons at energies of up to a few TeV to freely escape the system. The 
electromagnetic cascade is initiated by IC photons having much higher energies, up to ~ 10 4 
TeV, for which the 77-optical depth is much larger. The seed charges are provided by 
annihilation of MeV photons from the RIAF. It is found that the gap width is not smaller 
than 0.01r s if the density of seed charges is below the Goldreich- Julian (GJ) value. The 
luminosity of the VHE photons produced in the gap can account for the TeV luminosity 
observed by H.E.S.S. Any intermittencies of the cascade formation process would naturally 
lead to variability of both, the magnetospheric TeV emission and the resultant force-free 
flow, as observed. A schematic illustration of the model is presented in figure [TJ 



The strength of the magnetic field in a black hole magnetosphere is limited by the rate 
at which matter is accreted into the black hole. At sufficiently low accretion rates the flow 
becomes radiative inefficient (RIAF) and the electron temperature in the inner region of 
the RIAF may exceed m e c 2 (Narayan & Yi 1995, hereafter NY95). Bremsstrahlung cooling 
then gives rise to emission of soft gamma-ray photons that annihilate in the magnetosphere, 
leading to injection of charges on open magnetic field lines. For sufficiently high annihilation 
rate the resultant charge density can exceed the GJ value everywhere, keeping the magne- 
tosphere force-free. At lower annihilation rates the magnetosphere will be starved. In the 
latter case, the charge density produced by the annihilating photons from the RIAF defines 
the multiplicity required to establish a force-free outflow above the vacuum gap. 

In what follows we give a crude estimate of the ratio of the charge density produced by 
the RIAF and the GJ density, assuming that the accreting gas is in equipartition (/3 = 0.5) 
with the magnetic field. 

Henceforth, the accretion rate is measured in units of the Eddington rate, m = M / M^dd-, 
where the Eddington accretion rate is defined as M E dd = L Edd /r]ffC 2 = 10 27 M 9 gr s _1 , with 
r)ff = 0.1 adopted. The ion temperature of the accreted gas is close to virial, reaching 
Tj ~ 10 12 K at r = r s . The ion density is given by 



2. 



The role of a radiative inefficient flow 
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adopting v r = 0.1c(r/r s ) l ^ 2 for the RIAF radial velocity. The Thompson optical depth is 

r ~ ltfm{r/r s )- l/2 . (2) 

Typically r << 1 for accretion rates m « m CT it, where m cr ; t is the critical rate below which 
RIAF can exist. 

We assume that the equipartition magnetic pressure is half the gas pressure, viz., 
B 2 /8ir = 0.5piC 2 , where c s ~ c/a/3 x (r/r s ) -1 / 2 is the sound speed. This yields 

B-4x 10 4 (m/M 9 ) 1/2 (r/r s )- 5/4 G. (3) 

The numerical values in the expressions for n« and B are in rough agreement with the results 
of NY95 for a viscosity parameter a = 0.3, adiabatic index of 4/3, and radiative efficiency 
1 — / « 1, where / is the advection parameter, as defined in NY95. 

At radii r < 10 3 r s the electron-ion coupling becomes weak by virtue of rapid cooling. 
In this region electron cooling is dominated by synchrotron emission from thermal electrons. 
For accretion rates near the critical rate m cr ; t the electron temperature saturates at T e ~ a 
few times 10 9 K. At such temperatures only photons emitted by electrons at the tail of the 
thermal distribution have energies above the pair production threshold. However, for much 
lower accretion rates, fn « rn cvltl the electron temperature approaches 10 10 K, and thermal 
photons can annihilate. 

To estimate the annihilation rate at highly sub-critical accretion rates, we adopt the 
cooling functions for electron- ion and electron-electron Bremsstrahlung, q ei and q ee , from 
NY95. The total cooling rate per unit volume is then given to a good approximation by 
Qff = <?ee + Qei — 7-5 x 10~ 22 n 2 9 e erg s -1 cm -3 at electron temperatures 9 e = kT e /m e c 2 £ 1. 
The numerical value corresponds to the choice 9 e = 1 in the logarithmic term in Eqs. (3.6) 
and (3.8) of NY95. Since the RIAF is optically thin, pair production does not affect the 
leptonic content of the gas, and one can safely assume n e = rii. By employing Equation (CQ) 
one arrives at 

q ff ~ 1.8 x 10 2 6 e m 2 M~ 2 (r/r s )~ 3 erg s" 1 cm~ 3 . (4) 

The free-free luminosity emitted by the RIAF is Lff = J qffd 3 r ~ 2irr 3 qff ln(r/r s ), from 
which we readily obtain the number density of MeV photons in the magnetosphere: 

= q ff 2^Hr/r s ) ^ 0.2 g ff r 3 ^ ^ x 1Qll . 2 , 
2ncr 2 e~f cr i e ~t 

where e 7 = 39 e (m e c 2 ). The production rate of e pairs inside the magnetosphere due to 
77-annihilation is approximately o" 77 n 2 / c(47r/3)r 3 . In steady state this rate is balanced by 
the escape rate, roughly 4irr 2 n±c. Equating the two rates one has 

n± = a 77 n 2 r s /3 - 3x 10 u m 4 M^ 1 cm" 3 . (6) 
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The GJ density can be related to the accretion rate rh through Eq. (j3]) with r = r s : 

n GJ = ^ = 5x 10- 2 m^ 2 M~ 3/2 cm- 3 . (7) 

Thus we obtain 

n±/n GJ ~ 6 x 10 12 m 7/2 M 9 1/2 . (8) 

As seen, below a certain accretion rate, roughly m^,2x 10 -4 , injection of charges by this 
mechanism cannot provide complete screening of the magnetosphere. The ratio ncj/n± 
then defines the multiplicity required to produce a force-free flow. At higher accretion rates 
a vacuum gap may not exist. However, as noted above, the electron temperature decreases 
with increasing rh, eventually dropping below m e c 2 . Detailed ADAF calculations in Kerr 
spacetime (e.g., Manmoto, 2000; Li et al. 2009) indicate a cutoff in the emitted spectrum 
at around 9 e ~ 1 even for fn « m C rit- For 9 e < 1 the number density of MeV photons is 
suppressed by a factor of roughly 6 exp(l/6 l e ), and n± oc n 2 by a factor of [6 exp(l/6 l e )] 2 . Thus, 
there is an uncertainty of two to three orders of magnitudes in our estimate of the pair density 
n±. This large uncertainty and the strong dependence of charge injection on accretion rate 
motivates detailed self-consistent calculations. A complete, self-consistent treatment should 
also account for general relativistic effects, that are important close enough to the horizon. 

The strong dependence on accretion rate suggests that a gap may form during periods 
of low accretion, so that emission from the gap may be intermittent. 



3. TeV emission and pair production from charges accelerating in a starved 

magnetosphere 

We consider a rapidly rotating black hole of mass M = 1O 9 M 9 M embedded in an 
ambient radiation field. The radiation source is characterized by a luminosity Ld = 10 41 Z/4i 
erg/s, a size Rd = lZr s = 3 x 10 14 M 9 7?. cm, and SED (i.e., vF v ) peak energy eo = leo eV. The 
energy density of this radiation field is u s ~ "iL^{Rr 2 M^ 2 erg cm -3 , and the corresponding 
number density of photons at the peak (measured over a logarithmic energy interval) is 

n s (e ) ~ 1.8 x 10 12 L 41 ^- 2 M 9 - 2 e 1 cm" 3 . (9) 



3.1. Curvature and inverse Compton emission 



The electric potential difference across a gap of height h generated by a maximally 
rotating black hole, can be expressed as 

AV = 1.7 x 10 21 B A M 9 (h/r s ) 2 Volts. (10) 
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Charges accelerating in the gap will quickly reach a terminal Lorentz factor at which radiative 
losses balance energy gain, viz., qAV = Pt = Ph/c (Levinson 2000), where the net energy 
loss rate ,P = P cur + Pic, is the sum of curvature losses 

p- - (id 

3 p 2 

here p is the curvature radius of magnetic field lines, and inverse Compton (Thomson) losses 

Pic = o TC1 2 u s . (12) 

Equating AV and P CU rv and using Eqs. ([TO]) and ffTTj) . one finds that curvature radiation 
limits the Lorentz factor of emitting electrons (positrons) to 

lcur = 5 x lO^X^V^V^) 172 - (13) 

If the Thomson regime applies, inverse Compton scattering of electrons on the ambient 
photons limits the Lorentz factor to 

lie = (j^) V2 = 2 x 10 9 5 4 1/2 M 9 L- 1/2 ^( h/r s y'\ (14) 



Comparing Eqs. (1131) and (114jl . it is seen that particle losses are dominated by inverse 
Compton scattering, viz., jjc < jcur, if 

L 41 /TZ 2 > 1.6 x 10- 3 J B 4 1/2 M 9 (Vr s ) 1/2 (p/r s )- 1 . (15) 
The spectrum of curvature emission peaks at an energy 

e c „ = < 10B¥ 4 MV\h/r s )V\p/r s y/ 2 TeV, (16) 

l p 

where 7 max = min(7 cur , jic)- The corresponding number of curvature photons emitted by a 
single particle is, to a good approximation, 

iV 7 = P cur h/(ce cr , max ) > 4 x 10 8 s1 /4 M 9 1/2 (V^) 5/4 (p/^) -1/2 (1 + fr\ (17) 

where / = Pic /Pcur denotes the ratio of IC and curvature loss rates. Likewise, the maximum 
energy of IC photons is 

eic, max = m e c 2 lmax < 10 3 J B] /2 M 9 L 4 - 1 1/2 ^(/i/r s ) 1 / 2 TeV. (18) 



-7- 



The number of IC photons emitted by a single electron depends on the spectrum of the 
target radiation field. A crude estimate gives 

iV 7 = P Ic h/{ce IC , max ) > 10 e - 2 Bl /2 Ll{ 2 n-\h/r s )V 2 (l + (19) 

but the actual number may be much larger, roughly by a factor of m e c 2 / (j max hi/ s o) if 
the spectral energy distribution of the target radiation field peaks at a frequency v s q that 
satisfies hu s0 « m e c 2 . Inside the gap the field aligned electric field, E\\, is unscreened. 
This implies that the charge density on magnetic field lines must not exceed the GJ value, 
no j = QB cos8/2nec. The maximum gamma-ray power that can be produced by particles 
accelerating in the gap, regardless of the specific gamma ray production mechanism, is thus 

L i = J n GJ c(eAV)2nr 2 d9 ~ 3 x 10 47 'r]BlM 2 (h/r s ) 2 erg s _1 , (20) 

where rj is a geometrical factor. It is worth noting that the ratio L^/Lbz, where Lbz is the 
maximum BZ power that can be extracted by a force-free flow, scales as (h/r s ) 2 . If the pair 
density is well below the GJ density, then the gamma-ray power in Eq. ( 120]) is reduced by a 
factor of n±/ricj- 



3.2. 77-opacity and pair multiplicity 

The soft photon energy that corresponds to the pair production threshold is e t h r = 
m 2 c A /e 1 , which for peak curvature photons is 

e thr = 2.5 x 10- 2 £ 4 - 3/ X~ 1/2 (V^r 3/ Vr s r 1/2 eV, (21) 
and for IC photons is 

ethr = ^ = 2.5x l^B^M^L^K-Hh/rsY 1 ' 2 eV. (22) 
lie 

Let us denote by C( e s) the dimensionless energy distribution of soft photons, specifically 
the ratio of the number density at some energy e s and the number density of peak photons: 

n s {e s ) = C(e s )n s (e ). (23) 

Then, invoking <r 77 = 0.2o"t and using Eq. (jUJ), the pair-production optical depth to infinity 
at energy e s can be expressed as 

r 77 = a 11 n s (e s )R d ~ 70((e s )L 41 TZ~ 1 M g 1 e \ e s > e thr . (24) 
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The optical depth across the gap is smaller by a factor h/Rd, where we make the reasonable 
assumption that h < Rd- Note that the optical depth depends implicitly on the gap size 
h/r s through the value of C( e thr) at the threshold energy. 

Given Rd, Ld, Mg and B4 one can employ the above conditions in a self-consistent 
manner to solve for h/r s and 7 max . Let dN 1 /de 1 denote the number of photons per unit 
energy emitted by a single particle in the energy interval (e 7 , e 7 + de 7 ). The gap multiplicity 
M. can be defined as 

M = I (1 - e HV^K)]) ^l de ^ (25) 

and the gap size is then determined from the condition M. = ncj/n±, where n± is the density 
of seed charges injected into the open magnetosphere (see Eq. [8j). If the latter condition 
cannot be satisfied it means that the gap must be fully restored, with h/r s ~ 1. An outflow 
may still be established above the gap if the size of the external radiation source is large 
enough, Rd » r s ~ h. In such a case the gamma ray photons emitted by the accelerating 
charges in the gap will interact with target photons at larger radii, r > h, ensuing pair 
cascades there. The total multiplicity can then be computed from Eq. (|25|) upon setting 
h/Rd = 1 in the exponent. A force free flow will form above the starved magnetosphere 
provided M. > n±/ricj- 



4. Application to M87 and Sgr A* 
4.1. M87 

For this source we adopt a black hole mass of Mg = 4. The assumption that the M87 
jet is extracted magnetically by a BZ mechanism imposes a lower limit on the magnetic 
field in the vicinity of the horizon, and on the accretion rate. Various estimates (see, e.g., a 
summary in Li et al. 2009) yield a range of Lj ~ (10 43 — 10 44 ) erg/s for the power of the M87 
jet. Equating Lj with the BZ power that can be extracted from a Kerr black hole having a 
mass Mg = 4 and specific angular momentum a, 

L BZ = — a 2 B 2 r 2 s c ~ 3 x 10 47 ea 2 m erg s~\ (26) 
64 

where e ^ 1 is an efficiency factor that depends on field topology and other details, implies an 
accretion rate fn > 10 _4 e _1 a -2 . For a maximally rotating hole, a = 1, the implied accretion 
rate is ~ 10~ 4 , so that the magnetic field strength close to r s becomes B ~ 200 G. For an 
ADAF model with Lf, ~ 10 m 2 LEdd (NY95), this value for the accretion rate is consistent 
with limits on the bolometric luminosity measured in M87, L& < 10 42 erg s -1 (Owen et al. 
2000). Note that the Bondi accretion rate in M87 inferred from Chandra measurements is 
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tub ~ 1.6 x 1(T 3 (Di Matteo et al. 2003). The radial inflow velocity in a RIAF is about a 
factor a lower than in a Bondi flow, where a ~ 0.1 is the viscosity parameter, so that one 
expects 7Ti ~ arfiB (Narayan 2002), which would be consistent with the estimate above. 

The condition (jSj) formally yields n±/ncj = 0.12 (m/10 -4 ) 7 / 2 . However, as explained 
above this estimate is highly uncertain by virtue of the sensitive dependence of the injected 
pair density n± on the temperature of the accreting gas, on the geometry of the emitting 
region, and GR effects. It seems likely, though, that a multiplicity of 10 2 to 10 3 should in 
any case be sufficient for complete screening. 

If sufficient charges cannot be provided by emission from the RIAF to ensure that 
n± > n GJ , then a gap will exist, where the seed charges produced by annihilation of the 
RIAF photons accelerate. These seed charges will emit TeV photons that can initiate pair 
cascades via interaction with the ambient radiation, as discussed above. However, observa- 
tions constrain the 77-opacity at photon energies of a few TeV (e.g., Aharonian et al. 2006), 
and so pair cascades can only be initiated by photons having energies well above 10 TeV. 
In what follows, we employ some of the results derived in the preceding section to impose 
constraints on the maximum multiplicity and on the gap height. 

Provided the RIAF SED is not significantly modified by the presence of an outflow, 
the emerging spectrum ranges from radio to X-ray energies. The radio to sub-mm regime 
is produced by synchrotron radiation from relativistic thermal electrons and dominated by 
emission from small radii. Compton up-scattering of the synchrotron photons fills in the 
sub-mm to hard X-ray regime, with most of the emission coming from the inner ( ^ 10 r s ) 
region (e.g., Manmoto et al. 1997). For high rh, the Compton component is roughly a 
power law, while for low m distinct Compton peaks appear. This is related to the increase 
in electron temperature with decreasing accretion rates. At low rh, the X-ray spectrum is 
dominated by bremsstrahlung emission with a non- negligible contribution from large radii. 

The synchrotron emission is self-absorbed, resulting in a black-body spectrum up to 
frequency u sa , where it becomes optically-thin. In terms of ADAF parameters, one has 

u sa (r) ~ 9 x 10 10 M 9 - 1/2 (m/10~ 3 a) 1/2 ^9 {r s /rf 4 Hz , (27) 

where T e 9 = T e /10 9 K is the electron temperature (NY95; Mahadevan 1997). The luminosity 
then becomes (NY95; Mahadevan 1997) 

u sa L Vsa c 9 x HPMf (^L- J u 7 JZ n T$* erg/s, (28) 

where z/ n = u/10 11 Hz. As T e varies only slowly with r in the inner region of an ADAF, the 
highest synchrotron emission frequency comes from the innermost part (r ~ r s ). For M87 
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with Mg = 4, rh ~ amg and T ei g ~ 5 (cf. Mahadevan 1997), the synchrotron peak is located 
in the sub-mm band at 

v P = v.sa(r s ) ~ 10 12 Hz (29) 
and carries a characteristic power of 

L p = u p L Up ~ 10 41 erg/s . (30) 

This output would still comply to the upper limit imposed by Spitzer (MIPS) observations 
(several arcsec resolution), indicating that the (nuclear region plus galaxy) flux in M87 at 
2 x 10 12 Hz is ~ (3 - 4) x 10 41 erg/s (Shi et al. 2007). 

Compton up-scattering by a thermal distribution of relativistic electrons of small scat- 
tering depth r ~ 10 _2 (m/10~ 4 ) (rg/r) 1 ' 2 , Eq. (|2J), is characterized by a mean amplification 
factor per scattering of A = 1 + 48 + 166 2 . For ADAF temperature ranges of interest, typ- 
ically 2 < A 10 2 . The frequency of the first (once-scattered) Compton peak is roughly 
v c \ ~ Au p , and the power at this frequency is Lei ~ r] g L p rA where r) g ~ 0.5 is a geometrical 
factor. 

Pair-production implies that TeV photons of energy e produced in the magnetosphere 
will interact most efficiently with ambient soft photons of energy e s ~ (1 TeV/e) eV. For 
moderate A ~ 10, the target soft photon field is dominated by the fraction of sub-mm photons 
scattered twice (2nd Compton peak) with relative peak flux L C 2/L P ~ A 2 r 2 ~ 10~ 2 . Eq. (T2"4"|) 
then suggests that r 77 (e) < 1 up to VHE photon energies of several TeV. 

If, on the other hand, A ^> 10, then the energy of the once-Compton-scattered syn- 
chrotron photons can approach the infrared regime, v c ~ 10 12 A Hz, yielding a non-negligible 
background for 77-interaction with the TeV (curvature or IC) photons. Taking L 41 (e s ) ~ 
L C i ~ 0.3, Eq. (|2~4|) with 1Z ~ 1 formally implies that the pair-production optical depth for 
photons of e ~ 1 TeV is T 77 (e) ~ 5. However, this estimate is again somewhat uncertain by 
virtue of its sensitive dependence on accretion rate: 

The electron temperature in the inner region of an ADAF is only weakly dependent 
on accretion rate and at low rh roughly follows T e oc rh~ q with < q ^ 0.2 (Manmoto 
et al. 1997; Mahadevan 1997). The flux in the synchrotron peak, on the other hand, is 
highly dependent on the electron temperature and approximately varies as L p oc m 3 / 2+3 / 4 T e 7 
(Mahadevan 1997). For high enough temperatures A oc T 2 . This suggests that the Compton 
power, and thus the optical depth may roughly vary with accretion rate as 

r 77 (TeV) oc Lei oc tL p A oc rh 2 . (31) 

An even more sensitive dependence r 77 (e) oc m 2 5-4 has been estimated based on detailed 
RIAF modelling of M87 (Li et al. 2009). Hence, only moderate changes in the accretion are 
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required to make the source transparent to pair absorption up to photon energies of a few 
TeV as expected from IACT observations (e.g., Aharonian et al. 2006; Acciari et al. 2009). 
On the other hand, for energies well above 10 TeV severe absorption seems unavoidable. 

The RIAF photons also provide a suitable target field for inverse Compton interac- 
tions with the energetic charged particles accelerated in the magnetosphere. As the energy 
loss rate for a single electron in the extreme Klein-Nishina limit is only weakly (logarith- 
mically) dependent on the electron Lorentz factor, IC interactions (Thomson regime) with 
the radio-sub-mm disk photons (with the peak photons coming from r ~ r s ) becomes most 
constraining. IC interaction with photons of energy e p = hv p occurs in the Thomson regime 
for electron Lorentz factors up to 7 e ~ 10 s . If the gap size is not too small (h/r s ^ 10~ 2 ), 
electrons can get accelerated beyond this point (Pic(le) h/c < eAV). Taking B A ~ 0.02 
and 71 ~ 1 for M87, IC scattering f[T4"|) would limit achievable electron Lorentz factors to 
7ic ~ 10 9 L t l( 2 (h/r s ) 1 ^ 2 , where L t ^i = L t /10 41 erg s _1 is the luminosity of the target field 
for which Compton scattering occurs in the Thomson regime ihv t < m e c 2 f^). If one takes 

first approximation 

L t (u t ) ~ L p f^j , (32) 
as inferred from Eq. (I2"8"j) . IC interactions would limit electrons Lorentz factors to 

7ic^4xl0 n L;jf(Vr s ) 5/3 . (33) 

Note, however, that as vL v oc z/ 7 / 5 , with v oc r~ 5 / 4 , the dominant contribution to L t comes 
from larger scales, TZ t = Tl.{v t ) > 1, making the constraint even weaker, roughly by a factor 
of Jl] ^ 3 . Hence, if the gap is not too small, the energy loss of an electron will be dominated 
by curvature emission ( 1131) . limiting achievable Lorentz factors to 

/ h \ 1/4 

7c-3.8x 10 10 f-J , (34) 

where a curvature radius p = r s has been invoked for the magnetic field lines. 

Equation (TIB"]) then indicates that the peak energy of the curvature spectrum is limited 
to ~ 50m 3//8 TeV, which is below 1 TeV for m £ 10 -4 . However, the maximum energy of 
IC photons up-scattered by the accelerated electrons exceeds 7 c m e c 2 ~ 10 4 TeV. Given the 
expected spectral index of the radio disk spectrum in M87, cf. Eq. fl28l) . the average IC power 
emitted by a single electron in the Thomson regime; that is, due to scattering of photons 
having v < v t = m e c 2 / , ~y c , is 

2 j / 2 \ 7/5 

-3xl0 3 L Pi41 ^ 2 (Vr s ) 3/2 ° erg/s. (35) 

2nR z d {v t ) \>y c his p J 
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These photons originate from a radius 1Z — IZt » 1- The IC power emitted in the KN 
regime (for target photons at u > u t ) is reduced by a factor of ln(i/ p /z/ t )(z/ t /z/ p ) 2 : 

P IC (v p ) ~ f c llf\ Mv p /is t )(is t /is p ) 2 ~ 6 x 10 2 L Pi41 ^; 2 erg/s, (36) 

and depends only weakly (logarithmically) on the gap height. These photons are emit- 
ted from 71 = 7l p ~ 1. The total power is Pic(jc) — Pic{ u t) + Piciy p ). The frac- 
tion of IC power to curvature power near the maximum energy, f m = Pic(jc) / Pcur(jc) — 
0.06L Pj41 ^- 2 (/i/r s )- 1 [l + 5{n p /TZ t ) 2 {h/r s ) 3 / 20 }, is likely to exceed a few percent. 

The optical depth for the up-scattered photons at the highest energy ejc — r y c m e c 2 ~ 10 4 
TeV, cf. Eq. §M§ with £(e s ) = (m e c 2 j^ c hv p ) 2 ^ , is likely to exceed r 77 (e 7C ) ~ 10 4 . Hence, 
these photons will be severely attenuated before escaping the source, producing pairs. If 
(h/r s ) is not too small, this will happen inside the gap. The arising multiplicity is of order 



N,(IC) = ~ 10 2 ' 9 %^ (-) [1 + 5(n p /n t ) 2 ( h/r s f 2 »} oc m 7 ^ , (37) 

using L p oc rh 9 ^~ 7q , 7 C oc m 1 ^ 8 . Typically, !Z t >> 1Z P , so that the second term in the square 
parentheses can be neglected. The condition that the photon number N-y(IC) should exceed 
the multiplicity required to screen out the magnetosphere, Ai = ncj/n±, where n± hereby 
refers to the seed charges, imposes a constraint on the gap height: 

h/r s >(^y 3 L;X 3 nr. (38) 

Note that the optical depth across the gap is ~ T 77 (e/c , )V Rd = T77 (ere) (V r s)^ _1 an d 
therefore larger than unity if h/r s ^ 10~ 4 1Z. 

For the above choice of parameters, equation fl20l) yields a gamma-ray luminosity of 
L 1 ~ 10 45 rj(h/r s ) 2 erg/s, which for r\ = 0.1, (h/r s ) 2 £ 10~ 3 is in excess of the observed 
TeV luminosity. The sensitive dependence of seed charges from the RIAF emission on the 
accretion rate strongly suggests that moderate changes in accretion rate will induce nonlinear 
variations of the gap height, and conceivably its occasional disappearance. This will result 
in intermittencies of the gap emission and the resultant outflow. 

While the curvature emission is anticipated to peak at sub-TeV energy (although any 
spread in the curvature radius of magnetic field lines will smear out the peak), the emitted 
spectrum at higher energies should depend on the details of the pair cascade process. The 
cascading gamma rays will be emitted from their local photosphere, at an energy that varies 
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with radius. As demonstrated in Blandfrod & Levinson (1995), inhomogeneities can lead to 
a very broad spectrum even in case of mono-energetic injection of electrons. Thus, detailed 
calculations of the cascade process are required to determine the shape of spectrum emitted 
from the vicinity of the gap. 

4.2. Sgr A* 

Estimates of the accretion rate that accommodate the extraordinary low bolometric 
luminosity of Sgr A*, L b < 10~ 8 L Edd , are in the range m ~ (0.1 — 5) x 10~ 5 , adopting a black 
hole mass M 9 = 4 x 1(T 3 (e.g., Narayan 2002; Nayakshin 2005; Cuadra et al. 2006). VHE 
observations of the Sgr A* region indicate a gamma-ray spectrum that can be described by 
a power-law extending beyond 10 TeV above threshold of 165 GeV, with an associated VHE 
luminosity of a few 10 35 erg/s (Aharonian et al. 2004). The source of this emission, whether 
associated with the putative BH or with a different object, e.g., the PWN G359.95-0.04, is 
unresolved yet. From Eq. (JSJ) we obtain n±/riGj — 10~ 6 (m/10~ 5 ) 7 / 2 . Equation (1151) implies 
that energy losses are dominated by curvature emission as long as K > 2 (m/lO^ 5 )- 1 / 8 ^ 2 , 
where L^q = 1 is roughly the bolometric luminosity observed in Sgr A*, which is dominated 
by emission at frequencies around the SED peak, v, p ~ 10 12 Hz. The potential drop across 
the gap is limited to AV ~ 10 18 (m/10~ 5 ) 1 / 2 V, and from Eq. f[T3"j) we deduce that the 
Lorentz factor of accelerating electrons is 7 ~ 2 x 10 9 (m/10~ 5 )- 1/s . Thus, the peak energy 
of curvature photons is e cr ~ 200(m/10~ 5 ) 3 / 8 GeV, and the maximum energy of IC scattered 
photons eic is not expected to exceed ^m e c 2 ~ 10 3 TeV. The pair production optical depth 
at this energy is r 77 (e/c) ~ 10 2 7^ _1 , so the IC scattered photons are likely to be converted 
inside the gap. Equations © and JTSj) with L 41 = 10" 5 yield A 7 ~ lO^m/lO- 5 ) 1 / 4 ^ 1 , 
from which we infer that multiplicity may be high enough to allow the pair density to 
approach the GJ density if m > 10 -4 or so. For lower accretion rates we can safely assume 
h ~ r s . The present estimates of fn imply that a force-free flow cannot be established by 
this mechanism. 

The gap will, nonetheless, emit gamma-rays with a total luminosity 

L 7 ~ 10 35 (m/10- 5 ) 9/2 r? erg/s. (39) 

This emission is anticipated to peak in the sub- TeV band, at ~ 200 GeV. Some fraction may 
be emitted also at much higher energies. 
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5. Summary 

Annihilation of MeV photons that are emitted by a radiative inefficient flow in the 
vicinity of a supermassive black hole leads to injection of charges into the black hole mag- 
netosphere. The density of the injected charges depends sensitively on the accretion rate. 
At sufficiently large accretion rates the density may exceed the GJ density, giving rise to a 
complete screening of the field aligned electric field, and the formation of a force-free outflow. 
At lower accretion rates complete screening cannot be accomplished and a gap forms. The 
seed charges injected into the gap are quickly accelerated along magnetic field lines by the 
large potential drop generated by the rotating hole, until reaching a Lorentz factor at which 
energy gain is balanced by curvature and inverse Compton losses. The curvature and IC 
spectra peak in the TeV band. The interaction of the TeV photons thereby produced with 
the ambient radiation field may give rise to initiation of pair cascades, with a multiplicity 
that depends on the accretion rate. 

The proposed scenario is applicable to variable TeV-emitting Galactic Nuclei that ac- 
crete at sufficiently low rates. Out of the nearby under-luminous TeV candidate sources 
(Sgr A*, Cen A, M87, NGC 1399), M87 and Sgr A* are particularly interesting^] 

In case of M87 we estimate (with a large uncertainty) that annihilation of MeV photons 
from the RIAF cannot provide complete screening for accretion rates that correspond to the 
inferred jet power, and to a fit of the observed SED by an ADAF model. However, for a fully 
restored gap a large multiplicity, in excess of 10 3 , is expected owing to absorption of >> 10 
TeV photons that are produced through inverse Compton emission of the charges accelerating 
in the gap. The pair production opacity decreases with decreasing gamma-ray energy, and 
becomes sufficiently small below 10 TeV to allow photons at these energies to freely escape 
the system. We propose that the variable TeV emission detected by HESS was emitted from 
a magnetospheric gap located at the base of the VLBA jet. We note that rapid variability is 
naturally accounted for by this mechanism, as even moderate changes in accretion rate will 
lead to nonlinear fluctuations of the gap potential and the resultant TeV emission, owing to 
the sensitive dependence of the density of injected charges on accretion rate. The observed 
gamma-ray spectrum should depend on the spectrum of soft (scattered) photons and the pair 
cascade process, and can, in fact, be quite broad. Moreover, any spread in the curvature 
radius of magnetic field lines in the gap will lead to a spread in the maximum energy of 
accelerated pairs, which in turn broaden the spectrum further. Detailed calculations are 
needed in order to determine the shape of the emitted gamma-ray spectrum. Additional 



2 NGC 1399 has not yet been detected at TeV energies, while Cen A may not accrete at sufficiently low 
rates. 
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TeV photons may be emitted by the jet on larger scales, but this component should be more 
steady. It can be useful to trace and disentangle the spectrum of the variable component 
alone. 

Our analysis indicates that pair production alone would lead to a very high-sigma flow, 
which doesn't seem to be supported by the observations. Other processes, specifically 
magnetic field annihilation of a striped wind configuration by a Ray leigh- Taylor instabil- 
ity (Lyubarsky, 2010), might give rise to additional pair creation during the acceleration 
phase of the outflow. This, however, would still imply asymptotic Lorentz factor well in 
excess of that inferred from observations (albeit on much larger scales). Mass loading of 
the outflow may be accomplished through entrainment of surrounding matter, e.g., due to 
a non-dipolar magnetic field (McKinney & Blandford, 2009), or some instability of the wall 
jet interface. The resulting structure may consist of a fast (high-sigma) flow, surrounded by 
a slower sheath (e.g., Levinson & Eichler 1993). The very inner magnetosphere may still 
support a vacuum gap if entrainment occurs on scales larger than several r g . The details of 
such putative structures remain to be explored. 

In case of Sgr A* we find that for current estimates of the accretion rate the density of 
pairs in the magnetosphere is likely to be below the GJ value. Gamma-ray emission from 
the gap is expected with a total luminosity that depends sensitively on accretion rate, Eq. 
( 139]) . and an SED peak in the sub- TeV band. 

We thank John Kirk, Yuri Lyubarsky and Jon McKinney for enlightening discussions. 
AL acknowledges support by an ISF grant for the Israeli Center for High Energy Astro- 
physics. F.R. acknowledges support by a LEA Fellowship. 
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Fig. 1. — Schematic representation of the magnetosphere structure: A vacuum gap of height 
h < r s accelerates particles (electrons or positrons) to high Lorentz factors. The gap is 
exposed to soft radiation emitted from a source of size Rd- Curvature emission and inverse 
Compton scattering of ambient radiation produce VHE photons with a spectrum extending 
up to 10 4 TeV. Photons having energies below a few TeV can escape freely to infinity. 
Interactions of IC photons having energies well above 10 TeV with the ambient radiation 
initiate pair cascades just above the gap, leading to a large multiplicity. A force free outflow is 
established just above the pair formation front, and appears as the VLBA jet. Intermittencies 
of the cascade process, induced by modest changes in accretion rate, give rise to the variability 
of the TeV emission observed by HESS, and the fluctuations of the resulted force-free outflow, 
as indicated by the morphological changes of the VLBA jet. 
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